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Abstract 

The shell structure underlying shape changes in neutron-rich nuclei between AT = 20 and N = 28 
has been investigated by a novel application of the transient field technique to measure the first- 
excited state g factors in '^^S and '^'^S produced as fast radioactive beams. Details of the new 
methodology are presented. In both ^^S and there is a fine balance between the proton and 
neutron contributions to the magnetic moments. Shell model calculations which describe the 
level schemes and quadrupole properties of these nuclei also give a satisfactory explanation of the 
g factors. In ^^S the g factor is extremely sensitive to the occupation of the neutron P3/2 orbit 
above the = 28 shell gap as occupation of this orbit strongly affects the proton configuration. 
The g factor of deformed does not resemble that of a conventional collective nucleus because 
spin contributions are more important than usual. 

PACS numbers: 21.10.Ky,21.60.Cs,27.30.+t,27.40.+z,25.70.De 
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FIG. 1: (Color online) Region of the nuclear chart indicating the primary beam isotopes ^'^Ar and 
^^Ca used to produce the ^^S and ^'^S secondary beams (red). The = 20,28 closed shells are 
shown outlined in black. Stable isotopes are shaded (blue). 



I. INTRODUCTION 

In a recent Letter [l we presented the first application of a high-velocity transient- 
field (HVTF) technique 3] to measure the g factors of excited states of neutron-rich nuclei 
produced as fast radioactive beams. Questions on the nature and origins of deformation 
between = 20 and A^ = 28 were addressed by measuring the g factors of the 2^ states 
in 15822 and 16S24. Fig. [H shows the relevant part of the nuclear chart and indicates the 
primary beams used to produce the isotopes of interest. In the present paper we give a more 
comprehensive description of the experiment and discuss the interpretation in greater detail. 
A more technically-oriented discussion of the new technique will be published elsewhere jsj. 

Figure 121 shows the 2i excitation energies and the B{E2] Of —>■ 2i) values for the sulfur 
isotopes. It is apparent from the reduction in the excitation energies of the 2f states and 
the increase in B{E2) values that the neutron rich isotopes between A^ = 20 and A^ = 28 
develop collective features. Furthermore, the A^ = 28 nucleus 16S28 does not have the high 
2f energy and small B{E2) like the closed-shell A^ = 20 isotope '^gS2o- The data in Fig. |21 
together with additional information on the low-excitation level structures imply that the 
even sulfur isotopes between A^ = 20 and A^ = 28 undergo a transition from spherical at 
16^20! to prolate deformed in ^^824, and 16^26) and that the A^ = 28 nucleus 16^28 appears to 



exhibit collectivity of a vibrational character p, M, IH IS| ■ However the evolution of collective 
features in these nuclei has underlying causes that remain unclear. Some have argued that 
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FIG. 2: Previously-measured -£'(2"'') and B{E2) values for the sulfur isotopes ^. 

a weakening of the = 28 shell gap is important j^, 1^, while others have opined that 
the effect of adding neutrons to the f'j/2 orbit is primarily to reduce the proton Si/2-d^/2 
gap and that a weakening of the = 28 shell gap is not needed to explain the observed 
collectivity near ^^S [l^. There have been several theoretical studies discussing the erosion 
of the = 28 shell closure and the onset of deformation (e.g. Refs. H, ^ and references 
therein). To help resolve questions on the nature and origins of deformation between = 20 
and A^ = 28, we have used a novel technique to measure the g factors of the 2^ states in 
16^22 ana i6^24- 

The remainder of the paper is arranged as follows: Section |n] describes the experimental 
procedures including radioactive beam production, targets and apparatus. The experimental 
results and details of the data analysis are given in Section ITTTl Our shell model calculations 
are presented and discussed in Section HVl 
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II. EXPERIMENTAL PROCEDURES 



A. Overview 

The experiments reported here use transient fields to study neutron-rich nuclei produced 
as radioactive beams by fast fragmentation. In our approach the nuclear states of interest 
are excited and aligned by intermediate-energy Coulomb excitation The nucleus is then 
subjected to the transient field in a higher velocity regime than has been used previously 
for moment measurements, which causes the nuclear spin to precess. Finally, the nuclear 
precession angle, to which the g factor is proportional, is observed via the perturbed 7-ray 
angular correlation measured using a multi-detector array. It is important to note that the 
transient-field technique has sensitivity to the sign of the g factor, which in itself can be a 
distinguishing characteristic of the state under study, since the signs of the spin contributions 
to the proton and the neutron g factors are opposite. 

The transient field (TF) is a velocity-dependent magnetic hyperfine interaction expe- 
rienced by the nucleus of a swift ion as it traverses a magnetized ferromagnetic material 



14l [l5 |. For light ions {Z ^ 16) the maximum TF strength is reached when v = Zvq, i.e. 
the ion velocity v matches the i^'-shell electron velocity Zvq [vo = c/137, Bohr velocity). 
Since the transient field arises from polarized electrons carried by the moving ion its strength 
falls off as the ion velocity exceeds Zvq and the ion becomes fully stripped; a transient-field 
interaction will not occur for fast radioactive beams with energies near 100 MeV/nucleon 
until most of that energy is removed. Thus slowing the fast fragment beams to velocities 
where the transient field can be effective is an essential feature of the measurement. 



B. Radioactive beam production and properties 

The experiment was conducted at the Coupled Cyclotron Facility of the National Super- 
conducting Cyclotron Laboratory at Michigan State University. Secondary beams of ^^S and 
^°S were produced from 140 MeV/nucleon primary beams of ^°Ar and ^'^Ca, respectively, 
directed onto a ~ 1 g/cm^ ^Be fragmentation target at the entrance of the A1900 fragment 
separator Q|. An acrylic wedge degrader 971 mg/cm^ thick and a 0.5% momentum slit at 
the dispersive image of the A1900 were employed. The wedge degrader allowed the produc- 
tion of highly pure beams and also reduced the secondary beam energy to ~ 40 MeV/nucleon. 
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Further details of the radioactive beams are given in TableHl The ^^S (^°S) measurement ran 
for 81 (68) hours. The 40 MeV/nucleon beams were made incident upon a target consisting 
of contiguous layers of Au and Fe, 355 mg/cm^ thick and 110 mg/cm^ thick, respectively. 



TABLE I: Production and properties of radioactive beams. 





Primary beam 






Secondary beam 




Ion 


Intensity 


Ion 


E 


Intensity 


Purity 




(pnA) 




(MeV) 


(PPS) 


(%) 


40 


25 


38s 


1547.5 


2 • 10^ 


> 99 


48Ca 


15 


40s 


1582.5 


2 • 10^ 


> 95 



C. Target design and test 

Prior to the transient-field measurement a test run was performed to ensure that the fast 
fragment beams at ~ 40 MeV/nucleon would be slowed to the appropriate velocity regime 
in the target layers and emerge from the target for downstream detection without severe 
energy straggling. A 1540 MeV ^^S beam was used for this purpose. 

Three targets were mounted, in turn, on a rotating target ladder and the energy dis- 
tribution of the emerging ions was measured with a 980 /im PIN detector placed 41 cm 
downstream of the target. The targets were: (i) the 355 mg/cm^ Au foil alone, (ii) the Au 
(355 mg/cm^) + Fe (100 mg/cm^), and (iii) the Au (355 mg/cm^) + Fe (110 mg/cm^) target, 
which was subsequently chosen for the gf-factor measurement. The angle between the target 
and the beam was varied to effectively increase the target thickness. The transmission of the 
beam through the target and the energy distribution of the transmitted ions were measured 
using the PIN detector. These measurements indicated that the Ziegler (1985) stopping 
powers |l3| for sulfur in Au and Fe targets were accurate to about ±5%, although they tend 
to underestimate the stopping in the Au layer of the target and overestimate the stopping 
in the Fe layer. 

The effect of the energy-width of the radioactive beam on the spectrum of the emergent 
ions is examined in Fig. El which compares the particle spectrum measured in the test 
run for the Au + Fe (100 mg/cm^) target with Monte Carlo simulations made with the 
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code GKINT_MCDBLS jl8[. In the simulations the energy spread of the beam was treated 
as a Gaussian distribution with its width specified by the Full Width at Half Maximum 
(FWHM) . In the case shown in Fig. |21 the sulfur fragments emerged with energies in the 
range from ~ 80 MeV to ~ 200 MeV. A similar post-target energy distribution was observed 
in the g factor measurements where both the beam energy and the thickness of the Fe 
layer were increased slightly. It is evident from Fig. El that most of this energy-spread 
stems from the energy width of the radioactive beam, which is determined largely by the 
momentum acceptance, Ap/p = 0.5% (FWHM), set at the dispersive image slits of the 
A1900 spectrometer. 



D. Apparatus 

Figure ID shows the experimental arrangement. The radioactive beams were delivered 
onto the Au + Fe target, which had dimensions 30 x 30 mm^. The target was held between 
the pole tips of a compact electromagnet that provided a magnetic field of 0.11 T, suffi- 
cient to fully magnetize the Fe layer. As the iron layer of the target is much thicker than 
those typically used in transient-field gf-factor measurements, the magnetization of a piece 
was measured using the Rutgers Magnetometer [19]. Fields of 0.062 T were sufficient to 
ensure saturation. To minimize possible systematic errors, the external magnetic field was 
automatically reversed every 600 s. 

Table im summarizes the properties of the 2^ states and the key aspects of the energy loss 
of the sulfur beams in the target, applicable for the g factor measurements. The high-Z Au 
target layer serves to enhance the Coulomb excitation yield and slow the projectiles to under 
800 MeV, while the thick iron layer results in a long interaction time with the transient field, 
maximizing the spin precession. The energies (-Ej), with which the sulfur ions enter the iron 
layer, were calculated taking into account the energy loss measurements for the Au layer 
alone, while the energies (-Eg), with which the ions emerge from the iron layer into vacuum, 
were determined from the Doppler shifts observed in the ^f-factor measurements. 

Projectiles scattered forward out of the target were detected with a 15.24 cm diameter 
plastic scintillator phoswich detector placed 79.2 cm downstream of the target position. The 
phoswich detector consisted of a 750 /im layer of fast BC-400 scintillator and a 5.08 cm layer 
of slow BC-444 scintillator. The maximum scattering angle, 5.5°, limits the distance of 
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FIG. 3: (Color online) Spectrum of 1540 MeV ^'^S beam particles after passing through 355 mg/cm^ 
of Au followed by 100 mg/cm^ of Fe. The upper panel shows simulations that examine the effect 
of the energy-width of the ^'^S beam, which is assumed to have a Gaussian distribution about 1540 
MeV with a specified Full Width at Half Maximum (FWHM). These simulated spectra all have the 
same number of events. The experimental spectrum in the lower panel agrees with the simulations 
which assume the beam energy has a Gaussian distribution with a FWHM of 1.4% (i.e. a standard 
deviation of 0.5%). The equivalent momentum width is Ap/p = 0.7%. 



TABLE H: Nuclear parameters and reaction kinematics. B{E2)] = S(£^2;0+ 2f). {Ei^e) and 
{vi^e) are the ion kinetic energies and velocities at the entrance and exit of the iron layer. The 
effective transient-field interaction time ieff is evaluated for ions that decay after leaving the target. 



Isotope 


E(2+) 


B{E2) t 


M) 


{E^) 


{Ee) 


{vi/ZvQ) 








(keV) 


(e^fm^) 


(ps) 


(MeV) 


(MeV) 






(ps) 




1292 


235(30) 


4.9 


762 


123 


1.75 


0.71 


2.98 


40s 


904 


334(36) 


21 


782 


145 


1.73 


0.75 


2.99 
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FIG. 4: (Color online) Overhead view of the SeGA detectors and side view of the target chamber 
showing the magnet, target and phoswich detector. Projectiles were detected in the downstream 
phoswich detector. 

closest approach to near the nuclear interaction radius in both the Au and Fe target layers. 
Positioning the particle detector downstream also lowers the exposure of the 7-ray detectors 
to the radioactive decay of the projectiles. 

While the sulfur fragments do not penetrate beyond the fast scintillator, the slow scin- 
tillator helps discriminate against more penetrating radiation such as light ions produced in 
the secondary target, or accompanying the beam, and (3 decays in the phoswich from the 
decay of the implanted radioactive beam. Along with the fast-slow particle identification 
information from the phoswich, the particle time-of-flight was also recorded with respect 
to the cyclotron RF. Triggers from the radioactive decay of the beam particles, which are 
much lower in energy than the beam particles, were minimized by raising the threshold of 
the phoswich discriminator. 

For the ^'^S run, a circular Pb mask was placed on the phoswich detector to block particles 
in the range 0° < 6* < 2.5° and hence lower the count rate by excluding those scattering 
angles where the Rutherford cross section is large but the Coulomb excitation cross section 
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FIG. 5: Relation between beam axis, magnetic field axis, and the coordinate systems used in the 
present experiment and its analysis. 

is small. The mask helped reduce random particle-7 coincidences and pileup events in the 
phoswich detector. No mask was used for the ""^S run since the particle rate was low enough 
for pileup to be negligible. 

To detect de-excitation 7 rays, the target chamber was surrounded by 14 HPGe detectors 
of the Segmented Germanium Array (SeGA) 2^. The SeGA detectors were positioned with 
the crystal centers 24.5 cm from the target position. Six pairs of detectors were fixed at 
symmetric angles (±^,0) = (29°, 90°), (40°, 131°), (60°, 61°), (139°, 46°), (147°, 143°), and 
(151°, 90°), where 9 is the polar angle with respect to the beam axis and is the azimuthal 
angle measured from the vertical direction, which coincides with the magnetic field axis. 



?o make a connection 
IJ, Il5| , the locations 



Figure El indicates the coordinate frame and definitions of the angles, 
with the notation used for conventional transient-field measurements 
of the pair of detectors at the spherical polar angles {9, 0) and {9, + 180°) are written as 
{±9, 0); thus if {+9, 0) = (|^|, 0), then (-^, 0) = (|^|, + 180°). Each ±9 pair is in a plane 
that passes through the center of the target. 

Two additional detectors were placed at (6^,0) = (90°, 112°) and (24°, 0°) to assist the 
measurement of the angular correlation. All 14 detectors were used to measure the 7-ray 
angular correlations concurrently with the precessions. Since the precession angles are small, 
the unperturbed angular correlation can be reconstructed by adding the data for the two 
directions of the applied magnetic field. 

The positions of the SeGA detectors with respect to the target position were measured to 
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an accuracy of 2 mm using a theodolite system. This information was used to find the actual 
crystal locations in conjunction with the SeGA crystal segment positions measured by Miller 



et al. 



2l|. The 32-fold segmentation of the detectors improves the position determination 



of the 7 ray from the entire crystal length (8 cm) to near the segment length (1 cm), which 
is needed for Doppler corrections due to the high velocity of the emitting nuclei. The angle 
of 7-ray emission was deduced from the position of the detector segment that registered the 
highest energy deposition. With this algorithm the position resolution is near (but does not 
reach) the segment length. 

The master trigger for the data acquisition was set to record particle-7 coincidences as 
well as down-scaled particle singles. Thirty three energy signals were recorded for each 7- 
ray detector, corresponding to the 32 segments and the central contact. To differentiate 
between the signals from the fast and slow scintillator, the pulse from the phoswich detector 
was charge-integrated over the whole signal and over the tail (slow) part of the signal in 
separate QDC channels. Time differences were recorded between the phoswich detector and 
the 7-ray detectors. To assist with particle identification, the time-of-flight spectrum was 
also recorded for particles striking the phoswich detector. Finally, each event included a tag 
that identified the direction of the external magnetic field. 



III. EXPERIMENTAL RESULTS AND ANALYSIS 



A. Particle and 7-ray spectra 

A particle identification plot obtained with the phoswich detector during the ^^S mea- 
surement is shown in Fig. IHl This spectrum was obtained by plotting the integrated charge 
for the slow (or tail) part of the phoswich signal versus the integrated charge for the whole 
signal. Most of the intensity corresponds to the sulfur projectiles. It has a small tail compo- 
nent and therefore occurs in a band located along the full signal axis. The upsloping lines 
of intensity that deviate from this heavy-ion band are due to low-Z particles that punch 
through into the slow scintillator. Such ions accompany the beam, probably originating at 
the acrylic Image-2 wedge in the A1900 spectrometer; they are not predominantly produced 
in the secondary target. 

The spectrum in Fig. [3 shows the time-of-flight between the cyclotron RF and the 
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Full Signal (channel) 

FIG. 6: (Color online) Phoswich particle identification plot recorded during the ^^S run. Pulse- 
shape discrimination was obtained by plotting the integrated charge for the slow (tail) signal versus 
the full charge-integrated signal. 




Phoswich Full Signal (channel) 



FIG. 7: (Color online) Time-of-flight (between the cyclotron RF and the phoswich detector) versus 
the phoswich energy recorded during the ^^S run. The labeled features are (a) light particles 
registered in the slow scintillator, (b) ^^S ions in coincidence with 7 rays, and (c) down-scaled 
events. See text for further details. 

phoswich trigger pulse versus the phoswich energy signal (integration of the full pulse shape). 
Three distinct features are evident in this spectrum. The upper feature, labeled (a), cor- 
responds to events which cause a large amount of slow scintillation in the phoswich; they 
are the light particles discussed already. The middle region, labeled (b), corresponds to 
phoswich events triggered by sulfur projectiles that are in coincidence with 7 rays, and the 
lower region, labeled (c), is made up of downscaled phoswich events. (The downscaler mod- 
ule introduces an additional delay that shifts this group of events away from the phoswich-7 
coincidences.) 



12 



Gamma-ray spectra gated on sulfur recoils were produced and corrected for random 
coincidences by subtracting spectra gated on the appropriate regions of the particle-7 time 
spectra. Both the particle identification and time-of-flight information from the phoswich 
were used to select the events of interest. Spectra were also created for each 7-ray detector 
and for each direction of the magnetic field (up/down). Examples of the random-subtracted 
spectra for the lab-frame are given in the upper panels of Fig. |H] From the measured 
Doppler shift of the deexcitation 7 rays in the laboratory frame, the average after-target ion 
velocities were determined to be 0.083 c for ^^S and 0.088 c for ''^S, (i.e. v/Zvq = 0.71 and 
0.75, respectively). The velocity distribution of the exiting ''^S ions was also measured by 
shifting the phoswich detector by ±15 cm from its normal position and observing the change 
in the flight times of the projectiles. These procedures firmly establish that the sulfur ions 
were slowed through the peak of the TF strength at Zvq into the region where it has been 
well characterized ]^ [3| • 

Doppler-corrected spectra were also produced using the angular information from the 
SeGA detector segments and the particle energy information from the phoswich detector on 
an event-by-event basis, which is essential because of the spread in particle velocities. The 
lower panels of Fig. |H1 show examples of the Doppler-corrected 7-ray spectra. The peak to 
background ratio seems to be better for the '^°S measurement because (i) there are more 
counts in the ^'^S measurement, which gives an overall higher count baseline, and (ii) the 
combination of the higher 7-ray energy and shorter level lifetime in ^^S result in a broader 
Doppler-corrected peak. 

Most of the excited nuclei decay in vacuum after leaving the target material, however for 
^'^S a significant number also decay at higher velocities, extending up to the secondary beam 
velocity, whilst still within the target. The long Doppler tail becomes almost indistinguish- 
able from background at extreme forward and backward angles, but is clear at 6 = 40°, 60° 
and 90°. It was established through Doppler Broadened Line Shape (DBLS) calculations 
based on a Monte Carlo approach that the angular correlations and nuclear precessions 
can be determined accurately by an analysis of the vacuum-flight peak alone, without the 
need to include the Doppler tail. Examples of the DBLS calculations and comparisons with 
the experimental spectra are shown in Fig. El From these spectra it can be seen that the 
proportion of the 7-ray peak that corresponds to decays within the target is very small, 
and only weakly dependent on the detection angle. Note that the present calculations of 
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FIG. 8: Random-subtracted lab frame (top) and Doppler-corrected (bottom) 7-ray spectra at 
9 = 40° for ^^S (left) and ^'^S (right). The sulfur and iron peaks are labeled. The broad feature to 
the right of the 1292 keV line is mainly its Doppler tail, due to decays within the target, although 
there is a contribution from a contaminant line as well. The line at 1238 keV is the 4^2 transition 
emitted by ^^Fe nuclei at rest. 



the lineshapes do not account for the possibility that a 7-ray event may be assigned to the 
wrong segment of the detector. Having estabhshed that the correct angular correlation and 
gf-factor results are obtained by analyzing only the 7-ray peak, and excluding the tail, no 
attempt was made to obtain a quantitative fit to the observed 7-ray line shapes. Further 
details of these calculations will be given below and elsewhere 18[. 

The 2"*" peak areas averaged 925 counts/detector per field direction for ^^S and 400 
counts/detector per field direction for ^°S, in each of the six angle pairs of SeGA detec- 
tors used for extracting the precessions. 



B. Angular correlations and precessions 

In the rest frame of the nucleus, the perturbed angular correlation is given by 

w{e,(f),±Ae) 



k,q 



{2k + 1] 



X y^(|Ae|,±7r/2)r^, 



(1) 



where 



are the spherical polar angles at which the 7-ray is detected and ±AB is the 



transient-field precession angle for the two directions of the magnetic field; = 0, 2, 4 and q 
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FIG. 9: (Color online) Random-subtracted Doppler-corrected 7-ray spectra from the ^^S measure- 
ment are shown on the left. The 1238 keV line emitted by ^^Fe nucleus at rest in the lab frame 
are shifted towards the ^^S peak by the Doppler correction as the detection angle approaches 90°. 
The right panels show simulated spectra used to evaluate the impact of the Doppler Broadened 
Line Shape on the angular correlation and precession measurements. The insets show the energy 
range from 1150 to 1550 keV. These simulations, which do not include background, indicate the 
relative contributions of decays in the target (broad peak) and decays in vacuum (sharp peak), 
which overlap. The experimental peak is broader than that in the simulations because a simplified 
model is used for the 7-ray detectors. 

takes integer values such that —k < q < k. As shown in Fig. El the beam direction defines 
6 = and the magnetic-field direction defines = 0. The unperturbed angular correlation 
has symmetry about the beam axis and reduces to the 0-independent form: 

W{e) = 1 + E akQkGkPkicose). (2) 

k=2,4: 

The Ofc coefficients, which depend on the particle detector geometry, the spins of the initial 
and final nuclear states and the multipolarity of the 7-ray transition, can be evaluated from 
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the theory of Coulomb excitation [225. The correction factors for the finite sohd angles of 
the 7-ray detectors, Qk, are near unity in the present work. The deorientation coefficients, 
Gk, account for the effect of hyperfine fields experienced by ions that recoil into vacuum 
carrying atomic electrons. 

These recoil in vacuum effects were evaluated based on measured charge-state fractions 
for sulfur ions emerging from iron foils with energies between 92 and 236 MeV H- The 
large free- ion hyperfine interactions of hydrogen- and lithium-like ions {2^, which quickly 
reach hard-core values for nuclear lifetimes of a picosecond or more, are dominant. To a good 
approximation, the deorientation coefficients can be expressed as Ga; = 1 — {Qh + QLi)bk, 
where the fraction of ions with one and three electrons is Qh+Qu and bk = k{k+l) / (2/+!)^ 
for nuclear spin I and k = 2,4. Using the measured charge-state fractions gives G2 = 0.90 
and G4 = 0.65. 

The angular correlations were calculated with the program GKINT [25 1. For each layer of 
the target, the code GKINT performs integrals over the solid angle of the particle detector 
and over the energy loss of the beam in the target. The average coefficients, for example, 
are evaluated as 



imax 

P 



271 

("fc) = 1 I I ak{Op,E[z]) 



where 



a ■ ttgt) Jo 

X ^{9p,E[z]) sinOpdOpdz, (3) 



(a • ttgt) = 27r / / —{9p,E[z]) sin OpdOpdz, (4) 

where ^{9p, E[z]) denotes the Coulomb excitation cross section as a function of the scat- 
tering angle 6p and the projectile energy E[z], which varies with the depth z through the 
target layer. The integrals are evaluated numerically using Simpson's rule. To evaluate 
other average quantities of interest such as the average energy of excitation, or the average 
transient-field precession, the quantity to be averaged replaces in an expression of the 
same form as Eq. (jS)). 

Comparisons between the experimental and theoretical angular correlations are made 
in Fig E3 To indicate the effect of the Lorentz transformation (see e.g. j^, 0|) the 
angular correlations are shown in both the laboratory frame and the projectile frame. Good 
agreement was found between the calculated 7-ray angular correlations and the data. 
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FIG. 10: Angular correlations for ^^S and '^'^S. The left panels show the angular correlation in the 
frame of the projectile nucleus, while the right panels show the same angular correlations in the lab 
frame. The Lorentz boost causes a shift in the effective detection angle and an angle-dependent 
change in the solid angle. Data are normalized to the calculated angular correlation. The difference 
in anisotropy stems from the alignment produced by Coulomb excitation, which depends on the 



ratio of E{2'^) to the beam velocity (i.e. the adiabacity parameter, ^, see 
therein.) 



and references 



To describe the procedure for the extraction of the nuclear precession angle A6 it is 
helpful to begin with the 4 detectors in the plane through the target that is perpendicular 
to the magnetic field direction, as illustrated in Fig. ^2 For these detectors A© = A^^ and 
;he data analysis is identical to that in a conventional transient-field (/-factor measurement 



i^HH 



29|. The magnetic field causes a rotation of the angular correlation pattern 
around the x-axis, with positive precession angles for a right-handed co-ordinate frame in 
the direction indicated. Thus to first order in AG we have 

w ]i=w{eT^<S))^wie)T^<s>^, (5) 

do 

where the arrows indicate the direction of the external magnetic field. The precession angle 
is then determined from the double ratio of counts in a pair of detectors for each field 
direction. Considering, for example, the pair of detectors labeled 1 and 2 in Fig. ^2 the 
relevant experimental double ratio p is 
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Since the detector efficiencies and running times cancel out, p and A0 are related by 



where 

^^^^ - WlB 

It is conventional to define the 'effect', 



(8) 



e=(l-p)/(l + p), (9) 

so that 

AG = t/S. (10) 

These expressions apply only for the 4 detectors in the (f) = 90° (horizontal) plane. We now 
generalize to include the detector pairs that are not in this plane. By expanding Eq. to 
first order in A0 it can be shown that 

w Ti= w{e, 0, tA0) ~ w{e) t ag sin (11) 

do 

It follows that data analysis for a pair of detectors placed at the spherical polar angles {6, (p) 
and {6, + 180°) can proceed exactly as in the familiar case, but with the definition 

S{e,(f)) = S {6) sin (j). (12) 

Finally, it remains to note that the effect of the Lorentz transformation must be taken 
into account by evaluating 5* at the appropriate angle in the rest frame of the nucleus that 
corresponds to the laboratory detection angle 

The precession results are summarized in Table UTTl The two pairs of detectors at ^ = 29° 
and 151°, which were located in the horizontal plane perpendicular to the external magnetic 
field direction, near the angle of maximum slope of the 7-ray angular correlation, are most 
sensitive to the nuclear precession. 

To cover the uncertainties in the recoil in vacuum corrections and the possibility that some 
nuclear interference might affect the Coulomb excitation process for the (small) fraction of 
collisions that approach the nuclear interaction radius, an uncertainty of ±10% has been 
assigned to the S values. 
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FIG. 11: Schematic view of the experimental arrangement from above showing only the four SeGA 
detectors perpendicular to the magnetic field axis. The direction of positive precession angles, AG, 
is indicated. 



TABLE III: Precession results. 





detectors 






38s 






40s 




Pair 


e° 


^1 




e (xlO^) 


S 


A© (mrad) 


e (xlO^) 


S 


A0 (mrad) 


1 


29 


270 


90 


-50(20) 


+0.79(8) 


-64(27) 


-17(35) 


+0.86(9) 


-20(41) 


2 


151 


270 


90 


+26(23) 


-0.80(8) 


-33(29) 


-4(35) 


-0.88(9) 


+5(40) 


3 


139 


226 


46 


-10(23) 


-0.54(5) 


+19(43) 


+7(28) 


-0.58(6) 


-12(48) 


4 


147 


323 


143 


-7(27) 


-0.48(5) 


+15(57) 


-24(31) 


-0.52(5) 


+46(60) 


5 


40 


311 


131 


-41(19) 


+0.50(5) 


-83(39) 


+26(32) 


+0.53(5) 


+49(60) 


6 


60 


241 


61 


-18(15) 


+0.26(3) 


-69(58) 


+1(28) 


+0.27(3) 


+4(103) 


averag 


;e: 










-43(15) 






+5(21) 



C. Transient-field calibration and g factor results 
An evaluation of 



AQ/g = i-fi^/h) J B,,dt (13) 

is required to extract the experimental g factors. For light ions {Z < 16) traversing iron and 
gadolinium hosts at high velocity, the dependence of the TF strength on the ion velocity, w, 
and atomic number, Z, can be parametrized 0, as 



Bu{v, Z) = AZ^(t;/Zt;o)'e-5K^-o)*^ (14) 
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FIG. 12: (Color online) Transient field parametrization (solid line) for high-velocity light ions, 
:"rom Ref. Ml. Data points are measured transient-field strengths reported in the literature. See 

inn 

1, 115|, |30|| for further details and references. 

3ased parametrization that takes 







15 



30| and references therein). A 



where Vo = c/137 is the Bohr velocity. This is a model- 

n 

account of the known physics of the transient field (see |2, 
fit to data for iron hosts yielded A = 1.82(5) T with P = 3 The experimental data and 
the adopted parametrization are shown in Fig. El The velocity range sampled in the present 
experiments is also indicated. Although the data are sparse in the high velocity region, there 
can be no dispute about the general trend, and that the maximum TF strength is reached 
when the ion velocity matches the i^-shell electron velocity, v = Zvq. Also shown in Fig. IT^ 
are two alternative parametrizations of the field strength in the region v > Zvq chosen to 
give an indication of the uncertainty in the transient-field calibration. Compared with the 
adopted parametrization, these extrapolations give values of AQ/g that differ by ~ ±12%. 

Calculations of AQ/g were performed using the code GKINT to take into account the 
incoming and exiting ion velocities, the energy- and angle-dependent Coulomb excitation 
cross sections in both target layers, the excited-state lifetimes, and the parametrization of 
the TF strength in Eq. (fT^ . The results and the g factors extracted are given in Table HVl 

The g factor results are not very sensitive to the somewhat uncertain behavior of the 
transient field at the highest velocities because (i) the ions spend least time interacting with 
the TF at high velocity and (ii) the TF strength near 2Zvq is very small. Furthermore, the 
positive g factor in ^^S and the essentially null effect for are both firm observations, in- 
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TABLE IV: Experimental results. The spin precession per unit g factor, {AQ/g)caic, is eval- 
uated for ions that decay after leaving the target. The experimental g factor is given by 
g = AGcxp/(Ae/5)caic- 



Isotope 


(AG/5)calc 


AOexp 


9 




(mrad) 


(mrad) 




38s 


-330(41) 


-43(15) 


+0.13(5) 




-339(42) 


+5(21) 


-0.01(6) 



dependent of the transient-field strength. The g factor of '^^S is almost 3 standard deviations 
from zero. 

The experimental uncertainties assigned to the g factors are dominated by the statistical 
errors in the 7-ray count ratios, with small contributions from the angular correlation (10%) 
and transient field calibration (12%) added in quadrature. 

D. Monte Carlo simulations 

In the analysis of the data presented above an analytical formalism has been used to 
evaluate the average angular correlation coefficients, the transient-field precession per unit 
g factor, and other quantities of interest. 

To obtain a more detailed insight into the experiment and data analysis procedures it is 
helpful to use a Monte Carlo approach to model the experiment. The code GKINT_MCDBLS 
was written to generate event data for particle-gamma coincidences that can be sorted and 
correlated with calculated quantities such as the transient-field precession. The effect of 
the energy spread on the incident beam on the particle spectrum beyond the target was 
discussed above, and examples of the calculated Doppler broadened 7-ray line shapes were 
presented in Fig. IHl Here we examine the behavior of the precession angle and the angular 
correlation coefficients 02 and function of the detected 7-ray energy. 

Aside from the statistical uncertainties, which can be minimized by increasing the number 
of events in the simulation, the Monte Carlo approach allows for a rigorous treatment of the 
decay-in-fiight and vacuum deorientation effects. It also properly includes the contributions 
to the 7-ray peak that is analyzed to extract the g factor, stemming from excitation in the 



21 



Au and Fe layers of the target, and which cannot be separated experimentally. 

Figure ^1 shows how the angular correlation coefficients, 02 and 04, and the transient 
field precession, AQ/g, vary as a function of the 7-ray line shape in the laboratory frame. 
A detector at 29° in the laboratory frame was chosen for this comparison. The bump in the 
line shape near 1540 keV (lower panel) is due to the change in stopping powers as the ^^S 
ions pass into the Fe layer. As can be seen in the middle panel, there is no transient-field 
precession in the highest part of the Doppler tail (above the bump). This part of the line 
shape corresponds to ^^S ions that are excited and then decay in the Au layer of the target, 
before reaching the iron layer. As the sulfur ions slow and decay within the iron layer the 
transient-field precession AQ/g increases, reaching its maximum values for those ions that 
do not decay until they leave the iron layer and emerge into vacuum. The upper part of 
Fig. ^1 shows the variation of both ak and akGk, where Gk is the vacuum deorientation 
coefficient. Clearly, the vacuum deorientation effect applies only to those ions that decay 
in vacuum, and hence has an effect only in the peak part of the line shape. Aside from 
the vacuum deorientation effect, the variation in the magnitude of the ak coefficients as a 
function of the Doppler shift has its origin in the alignment produced by intermediate energy 
Coulomb excitation 2^, 2^. There is a pronounced dependence on the beam energy, such 
that the alignment is reduced for those projectiles that are excited after losing energy within 
the target. 

At present the variation of the 7-ray intensity across the acceptance of the detector 
is not included, and the response of the segmented Ge detector is treated approximately. 
For example, it is assumed that the hit is always assigned to the correct segment. These 
approximations are not significant in the present context. 

Because the ^^S nuclei can be excited in either the Au or Fe layers of the target, and many 
of them decay before they leave the target and emerge into vacuum, the analytic (GKINT) 
calculations of 026*2, a^G^ and AQ/g for the vacuum fiight peak required some simplifying 
approximations (mainly to estimate the relative contributions of excitation in the two tar- 
get layers). The Monte Carlo (GKINT_MCDBLS) calculations, however, enable a simple 
and rigorous evaluation of these quantities of interest. The disadvantage is that the Monte 
Carlo calculation is very time consuming. It was found that the difference between the two 
approaches was small in our case. As well as providing a deeper insight into the experiment, 
the Monte Carlo calculations therefore support our use of the approximate analytic calcu- 
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FIG. 13: (Color online) Monte Carlo simulations (10^ events). Lower panel: Doppler broadened 
7-ray line shape at 29° in the laboratory frame. An equivalent spectrum in the projectile frame is 
shown in the uppermost panel of Fig. |^ Much of the energy width in the laboratory frame stems 
from the spread in the projectile velocity, which must be Doppler-corrected event-by-event using the 
energy from the phoswich detector (section llIIA|) . Middle panel: The transient-field precession per 
unit g factor as a function of the emitted 7-ray energy. Upper panel: Angular correlation coefficients 
as a function of the emitted 7-ray energy. The effect of vacuum deorientation is indicated by the 
(red) line labeled UkGk- 



23 



lations, which introduce neghgible error compared with the statistical uncertainties in the 
experiment. 



IV. SHELL MODEL AND DISCUSSION 



Shell model calculations 



Shell model calculations were performed for ? fiS?n , 16S22 and igS24, and their isotones 
fgAr2o, i8Ar22 and fgAr24, using the code OXBASH and the sd-pf model space where (for 

> 20) valence protons are restricted to the sd shell and valence neutrons are restricted to 
the pf shell. The Hamiltonian was that developed in Ref. [s^ for neutron-rich nuclei around 

= 28, i.e. the SDPF-NR interaction [s^. These calculations reproduce the energies of 
the low-excitation states to within 200 keV. With standard effective charges of Cp ~ 1.5 
and Cn ~ 0.5 they also reproduce the measured B{E2) values. The relevant results for 
the 2i states are summarized in Table I3 where they are labeled SDPF. For the purposes 
of the following discussion, the B{E2) and Q{2~^) values are also presented in terms of 
the equivalent deformation parameter p. These electric quadrupole properties have been 

n 

calculated and discussed by Retamosa et al. [^M]- Although the shell model Hamiltonian 
has been improved since their work, the theoretical B{E2) and deformation values agree 
closely with the values in Table \V\ and their discussion remains relevant. 

The g factors of the 2f states were evaluated using the bare nucleon g factors. The 
calculated g factors are compared with experimental results in Fig. and Table IVH which 
also shows the orbital and spin contributions to the g factor originating from both protons 
and neutrons. As will be discussed below, the overall level of agreement between theory 
and experiment is satisfactory given the extreme sensitivity to configuration mixing and the 
near cancelation of proton and neutron contributions in the A^ = 22, 24 isotones. 

The main partitions of the 2^-state wavefunctions in ^^S and are indicated in Ta- 
ble IVIII Table lyTTTl shows the decomposition of the g{2~^) values into the contributions from 
each combination of single-particle orbitals for ^s.^og g^j^^ 38,40^^^, rjj^g off-diagonal contribu- 
tions from spin-orbit partner orbits can play an important role. In ^^S and they tend to 
quench the (diagonal) moment of the dominant configuration. 

Many authors have argued on experimental and theoretical grounds that ''^S (A^ = 24) is 
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TABLE V: Theoretical and experimental excitation energies, quadrupole moments, and B{E2) 



values. Experimental data are from 



Nuclide 


Model 


(MeV) 


(MeV) 


(e fm^) 


B{E2) 
(e^ fm^) 


B{E2) t<=^P 
(e2 fm^) 










SDPF 


3.426 


3.291 


-11.6 


141 


104(28) 


+0.21 


0.20 


0.17(2) 


38s 


SDPF 


1.531 


1.292 


-9.7 


268 


235(30) 


+0.17 


0.26 


0.25(2) 




SDF 


1.286 




+7.5 


158 




-0.13 


0.20 




40 s 


SDPF 


0.980 


0.904 


-19.3 


473 


334(36) 


+0.33 


0.34 


0.28(2) 




SDF 


1.052 




-9.8 


328 




+0.17 


0.28 




38 Ar 


SDPF 


2.018 


2.167 


+4.3 


178 


130(10) 


-0.07 


0.19 


0.16(1) 


40 


SDPF 


1.371 


1.461 


+8.5 


263 


330(40) 


-0.13 


0.22 


0.25(2) 




SDF 


0.966 




+12.8 


228 




-0.20 


0.21 




42 Ar 


SDPF 


1.243 


1.208 


+5.2 


360 


430(10) 


-0.08 


0.25 


0.28(3) 




SDF 


0.918 




+13.5 


317 




-0.20 


0.24 





deformed, some linking it to a weakening of the = 28 shell gap (see Refs. 



BBBay 



mm 

and references therein). To explore the role of excitations across the N = 28 shell 
gap, a further set of calculations was performed in which the neutrons were confined to the 
/7/2 orbit. The results of these calculations, labeled SDF, are given in Tables IVl and IVIl 



B. Discussion of results 



In the N = 20 isotones, ^^g and ^^Ar, the 2^ state is a pure proton excitation for 
our model space. In ^^S the first-excited state at 3.29 MeV is dominated (87%) by the 
proton (si/2C?3/2)2+ configuration for which g = +1.459. The 2f configuration in ^^Ar is 
predominantly (93%) TT{d3/2)l+ for which g = +0.083. This case demonstrates the extreme 
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TABLE VI: Theoretical g factors compared with experiment. 



Nuchde 


Model 




^proton 




^neutron 




^exp 






orbital 


spin 


total 








36 c 
16^20 


SDPF 


0.967 


0.187 


+1.154 





+1.154 




38 c 
16^22 


SDPF 
SDF 


0.225 
0.087 


0.073 
0.049 


+0.298 
+0.136 


-0.301 
-0.494 


-0.003 
-0.358 


+0.13(5) 


40 c 
16^24 


SDPF 
SDF 


0.225 
0.249 


0.051 
0.070 


+0.276 
+0.318 


-0.241 
-0.404 


+0.035 
-0.085 


-0.01(6) 


liArso 


SDPF 


1.151 


-0.842 


+0.308 





+0.308 


+0.24(12) ^ 


t^Ar22 


SDPF 
SDF 


0.311 
0.228 


-0.147 
-0.092 


+0.164 
+0.136 


-0.364 
-0.431 


-0.200 
-0.295 


-0.02(4) ^ 


fiAr24 


SDPF 
SDF 


0.263 
0.242 


-0.043 
-0.025 


+0.220 
+0.217 


-0.280 
-0.417 


-0.060 
-0.200 





"Presept|work. 
^Ref. 



m 



sensitivity of the magnetic moments to configuration mixing - the remaining 7% of the 
wavefunction raises the calculated moment to +0.308, in agreement with the experimental 
value js^. It also shows the importance of off-diagonal terms (see Table IVTTT|) as the diagonal 
contributions associated with the two main configurations account for only 60% of the total 
theoretical g factor. 

Two neutrons have been added to the pf shell in the N = 22 isotones ^'^S and ^°Ar. Since 
^^S is almost doubly magic, the initial expectation might be that the first-excited state of 
^^S would be dominated by the neutron fY/2 configuration weakly coupled to the '^^S core, 
resulting in a g factor near —0.3. An example of this type of behavior is the = 52 isotope 
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TABLE VII: Partitions > 5% contributing to the wavefunctions of the 2f states of and '^'^S in 
the SDPF model. 



nuchde 




proton orbit 




neutron orbit 




occupation (%) 




lc?5/2 


14/2 


2si/2 


1/7/2 


1/5/2 


2p3/2 


2pi/2 




38s 


6 


u 


2 


2 











26.89 




6 


-1 

1 


-1 

1 


2 











18.36 




6 


2 





2 











12.10 




6 


U 


2 


1 





-I 

1 





12.01 




6 


1 


1 


1 





1 





6.36 


40s 


6 


2 





4 











17.26 




6 


1 


1 


4 











16.22 




6 





2 


4 











9.30 




6 


2 





3 





1 


u 


y.io 




6 


1 


1 


3 





1 





8.37 




6 





2 


3 





1 





6.05 


TABLE VIII: Contributions to the g 


factors of the 2f states in ^^S, 


40 s 


Ar and ^°Ar for the 


SDPF model. 
















orbits 




38s 




40s 




38 


40 


TT 1^5/2 - 


lc^5/2 




0.046 




0.067 




0.044 


0.033 


TT 1(^5/2 - 


lt^3/2 




-0.015 




-0.003 




0.124 


0.003 


TT 1^3/2 - 


lf^3/2 




0.014 




0.014 




0.080 


0.021 


TT 2si/2 - 


2si/2 




0.254 




0.199 




0.061 


0.107 


1/7/2 - 


1/7/2 




-0.430 




-0.447 






-0.394 


2^ 1/7/2 - 


1/5/2 




0.039 




0.118 






-0.002 


2^ 1/5/2 - 


1/5/2 




0.009 




0.011 






0.005 


I' 2p3/2 - 


2p3/2 




0.054 




0.092 






0.018 


2p3/2 - 


2Pl/2 




0.026 




-0.016 






0.008 


U 2pi/2 - 


2Pl/2 




0.002 




0.000 






0.001 
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FIG. 14: (Color online) Theoretical g factors compared with experiment. The dashed line shows 
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the present shell model calculations. The previous results for '^^Ar and ^'^Ar are from Refs. 



and 



3q |. respectively. 



40^^52 



In contrast with 



37l |: the effect is also evident, but less pronounced, in 42M052 
this weak-coupling scenario, the near zero theoretical g factor and the small but positive 
experimental g factor of ^^S require additional proton excitations, which indicate strong 
coupling between protons and neutrons. It is relevant to note that strong proton-neutron 
coupling is considered one of the prerequisites for the onset of deformation. The effect on 
the proton configuration due to coupling with neutrons in orbits above the = 28 shell gap 
will be discussed below. 

For A^ = 22, 24 the shell model predicts a cancelation of the proton and neutron contribu- 
tions to the moment. Under these conditions, the description of the g factors is satisfactory. 
The dependence of the g{2~^) in ''^Ar on the basis space, the interaction, and the choice of 
effective nucleon g factors, has been investigated in Ref. j^. To estimate the impact of the 
use of effective g factors, we recalculated the g{2^) values in Table IVll adopting the effective 
gi and gs values employed by Stefanova et al. The effect was to make the 2"'"-state g factors 
more positive, by ~ 0.05 for the S isotopes and by ~ 0.1 for the Ar isotopes; the discrepancy 
between theory and experiment for ^^S and ''^Ar is therefore reduced, but not eliminated by 
the use of effective g factors. 

A comparison may also be made between ^^S and ''^S, and their isotones "^^Ca and ^^Ca, 
where the experimental g factors, 5'(2|^ -1-0.04(6) and g{2'^) = -1-0.16(3), are also far 
from that of the configuration 



ntly a 

significant collective component in the wavefunction that cannot readily be included in shell 
model calculations jssl Isol. Some remnant of this collective excitation may occur in ^^S and 
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^°S, although it is very much smaller than in the calcium isotopes. 

We conclude that the dominant components of the shell model wavefunctions are correct 
and give a satisfactory description of the measured g factors. As such, the shell model 
calculations provide a microscopic explanation of the development of quadrupole collectivity 
in these nuclei, which will now be discussed in more detail. 

The existence of deformation in nuclei has long been associated with strong interactions 
between a significant number of valence protons and neutrons, particularly in nuclei near 
the middle of a major shell. Without exception the deformed nuclei studied to date have 
g factors near the hydrodynamical limit, Z/ A, reflecting the strong coupling between protons 
and neutrons, and a magnetic moment dominated by the orbital motion of the proton charge 
with small contributions from the intrinsic magnetic moments of either the protons or the 
neutrons. Examples include ^'^Mg in the sd shell [4^, and ^°Cr in the pf shell 4l|. Robinson 



et al. 



42| have recently used the shell model to examine the relation between the B{E2), 



Q{2i) and the orbital magnetic dipole strength (scissors mode) in these nuclei. 

As noted above. Coulomb-excitation studies and the level scheme of suggest that it 
is deformed. Supporting this interpretation, the shell model calculations (Table |^ pre- 
dict consistent intrinsic quadrupole moments when derived from either the B{E2) or the 
quadrupole monrant, 0(2f), implying a prolate detojinatiou of /3 « +0.3, m agreement 
with the value deduced from the experimental B{E2) [5|. The near zero magnetic moment, 
however, does not conform to the usual collective model expectation of (7 ~ Z/A. Since the 
shell model calculations reproduce both the electric and magnetic properties of the 2i state 
they give insight into the reasons for this unprecedented magnetic behavior in an apparently 
deformed nucleus. 

The essential difference between the deformed neutron-rich sulfur isotopes and the de- 
formed nuclei previously encountered (i.e. either light nuclei with N = Z or heavier deformed 
nuclei) is that the spin contributions to the magnetic moments are relatively more impor- 
tant, especially for the neutrons. It can be seen from Table IVTl that the proton contributions 
to the g factors are dominated by the orbital component, as is usually the case for deformed 
nuclei, but the substantial neutron contributions originate entirely with the intrinsic spin. 
The main partitions of the shell- model wavefunctions summarized in Table IVTTl indicate that 
both ^®S and ^"^S have a dominant occupation of the neutron /7/2 orbit, for which g = —0.547 
(unless coupled to spin zero). The net neutron contribution to the g factor in '^'^S ("^^S) is 
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FIG. 15: Dominant proton partitions in S comparing calculations in which neutrons occupy 
the full pf shell (SDPF) with calculations in which they are restricted to the /7/2 shell (SDF). 
The g factors which represent the diagonal contributions of the most important configurations are 
indicated. 
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FIG. 16: As for Fig. [HI but for ^^S. 

therefore ~ 55% (~ 44%) of that of a pure f-jj2 neutron configuration. Those configurations 
with neutron excitations into the P3/2/5/2P1/2 shell mainly have a single occupation of the 
P3/2 orbit. Although these excitations tend to reduce the magnetic moment of the neutrons 
from the value of the pure configuration, they do not cancel the neutron magnetic mo- 
ment altogether because the g factor of the z/(/7/2 ® ^3/2)7 configuration coupled to spin 
J = 2 is —0.183, and the value increases in magnitude for other possible values of J . The 
off-diagonal terms in the Ml operator also quench the /7/2 neutron contribution to g{'2^)i 
as indicated in Table IVTTn 
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The role of neutron excitations across the = 28 shell gap can be examined further 
by comparing the calculations using the full SDPF model space with the truncated SDF 
calculations in which the neutrons were restricted to the /7/2 shell. Looking first at the E2 
properties presented in Table IVl it is apparent from the (5(2+) values that excitations into 
the P3/2/5/2P1/2 shell are needed to produce significant prolate deformations in ^^S and ^°S. 
Furthermore, in ^°S, the quadrupole moment and B{E2) are not consistent with the same 
intrinsic quadrupole deformation unless the neutrons can occupy the P3/2/5/2P1/2 shell. The 
wavefunctions in Table ETT] indicate that the development of quadrupole collectivity depends 
most sensitively upon the occupation of the 1^^3/2 orbit. 

Turning to the g factors, it is apparent that two effects are important. First, in all cases 
the neutron contribution to g{2~^) is quenched significantly when the P3/2/5/2P1/2 orbits are 
occupied. Second, the case of ^'^S shows a dramatic change in the g factor due to a relatively 
small neutron occupation of the ^3/2 orbit, highlighting that there is a strong coupling 
between the protons and neutrons. Specifically, a small occupation of the neutron P3/2 orbit, 
moves the theoretical g factor from —0.358 to —0.003, falling short of the experimental 
value by a small margin compared with the distance traveled. Fig. ^1 shows how the 
proton partition of the wavefunction for '^'^S changes significantly with the occupation of 
the P3/2/5/2P1/2 shell; there is a significant increase in the contribution of the proton 51/2^3/2 
configuration, which has a large g factor, and which has been proposed to drive deformation. 
In contrast. Fig. IT6l shows that for ^°S this proton configuration contributes strongly to the 
wavefunction whether excitations across the = 28 gap are allowed or not. Apparently the 
effect of a single neutron occupying the ps/2 orbit is diluted when three remain in the /7/2 
shell. 

Since the hydrodynamical collective models fail to explain the g factor of ^°S, the observed 
quadrupole collectivity in this region is better interpreted in terms of the symmetries of the 
shell model Hamiltonian. More specifically, the development of quadrupole collectivity can 
be hnked to the quasi-S'?7(3) symmetry identified by Zuker et al. [4^ and considered for 
the sulfur isotopes by Retamosa et al. 34]. In the neutron space the Aj = 2 orbits /7/2-P3/2 
develop the quasi-Sf/ (3) symmetry. As shown by the comparisons in Table |3 occupation 
of the P3/2 orbit is clearly essential if quadrupole collectivity is to emerge. Turning to the 
proton space, for the S isotopes of interest the proton ^5/2 orbit is essentially closed and only 
the and si/2 orbits are active (see Table IVTT|) . Although quasi-S'f/(3) cannot develop 
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for these two orbits, an approximate pseudo-5'[/(3) geometry does develop. When neutrons 
begin to fill the pf shell, the effective energies of the two proton orbits become increasingly 
degenerate, as is manifested by the narrowing of the gap between the lowest 3/2+ and 
1/2+ states in the odd-A isotopes of K, CI and P 10|, |4^. In the limit of degenerate ^3/2 
and Si/2 orbits the valence proton space has the geometry of pseudo-S't/(3). Within this 
framework, quadrupole collectivity develops in the neutron-rich sulfur isotopes because the 
proton number is optimal for quadrupole coherence, despite the fact that the ^3/2 - Si/2 
degeneracy is not reached 3^ . 

The shift in the effective single particle energies of the proton (^3/2 and si/2 orbits is 
strongly linked to the effect of the monopole component of the tensor term in the nucleon- 



nucleon interaction. Otsuka et al. 45| have shown that the effect of the monopole interaction 
between the proton ^3/2 and neutron /7/2 orbits is attractive, which narrows the gap between 
the proton (^3/2 and S1/2 states as more and more neutrons are added to the /7/2 shell. This 
effect of the monopole interaction is therefore more pronounced in than in ^^S, and 
may explain why the contribution of the proton partition 'i5/2'5i/2^3/2 ^^^^ sensitive to the 
neutron occupation of the P3/2 orbit in than in ^'^S. 

V. SUMMARY AND CONCLUSION 

In summary, we have developed a high-velocity transient-field technique to measure the 
g factors of excited states of neutron-rich nuclei produced as fast radioactive beams. The 
g factors of the first-excited states in the neutron-rich isotopes ^^S and have been mea- 
sured with sufficient precision to test and challenge shell model calculations. 

Keys to the success of these measurements on low intensity radioactive beams (10^ to 10^ 
pps) include (i) the use of intermediate energy Coulomb-excitation to align the nuclear spin, 
and (ii) exploitation of the high velocities of the fragment beams to maximize the transient- 
field precession angle. In these measurements the precession an gle p er unit q factor is ~ 10 
times that in conventional measurements on neighboring nuclei |35l {, 
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The nature and origin of the quadrupole collectivity that develops between = 20 and 
= 28, including the role of excitations across the = 28 shell gap, has been explored by 
comparing the measured g factors with shell model calculations. 

The case of ^^S highlights the role of strong proton-neutron interactions, showing that 
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neutron excitations across the N = 28 shell gap cause a large change in the proton configu- 
ration. Taken together with the results for ^°S, it seems that excitations across the = 28 
shell gap and the reduction in the proton Si/2-c?3/2 gap due to increasing neutron occupation 
of the fj/2 orbit both contribute to the development of collectivity. 

Since the hydrodynamical collective models fail to explain the g factor of ^°S, the observed 
quadrupole collectivity in this region should be interpreted in terms of the symmetries of 
the shell model Hamiltonian. The unusual combination of magnetic and electric properties 
in apparently comes about because relatively few particles are involved in the collective 
motion. These features reemphasize the unique mesoscopic nature of the nucleus. 

Looking to future applications of the high-velocity transient-field technique, it may be 
noted that the 2f state in the N = 20 nucleus '^2Mg2o has a similar excitation energy, lifetime 
and B{E2) to '^'^S. However the g factor in ^^Mg might be closer to that of a conventional 
collective nucleus since the = 20 shell closure is known to vanish far from stability. Indeed 
a .hen .ode. calculation tl. i,.e.act,o,. of Wa.b..to„, Beoke. a.d B.ow. U give. 
9 = 9p + 9n = 0.52 — 0.11 = +0.41, close to the collective value of Z/A = 0.375. 
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